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A CFD validation workshop for synthetic jets and turbulent separation control (CFD- 
VAL2004) was held in Williamsburg, Virginia in March 2004. Three cases were inves- 
tigated: synthetic jet into quiescent air, synthetic jet into a turbulent boundary layer 
crossflow, and flow over a hump model with no-flow-control, steady suction, and oscilla- 
tory control. This paper is a summary of the CFD results from the workshop. Although 
some detailed results are shown, mostly a broad viewpoint is taken, and the CFD state- 
of-the-art for predicting these types of flows is evaluated from a general point of view. 
Overall, for synthetic jets, CFD can only qualitatively predict the flow physics, but there 
is some uncertainty regarding how to best model the unsteady boundary conditions from 
the experiment consistently. As a result, there is wide variation among CFD results. For 
the hump flow, CFD as a whole is capable of predicting many of the particulars of this 
flow provided that tunnel blockage is accounted for, but the length of the separated region 
compared to experimental results is consistently overpredicted. 
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I. Introduction 


The study of separation control through the use of blowing, suction, and synthetic jets has been 
an active field of research for some time now. 1 Many experimental and theoretical papers (only 
a few of which are referenced here) have been published on the subject. 2-6 The review article by 
Greenblatt and Wygnanski 7 lists many additional references. There has also been a recent increase 
of attention paid to the field by the American Institute of Aeronautics and Astronautics (AIAA), 
which held its first Flow Control Conference in June 2002 in St. Louis, Missouri, and now sponsors 
it as a continuing biennial series. 

Computational methods have also been used extensively for computing synthetic jet flows. 
Only a few examples are referenced here. 8-13 However, to our knowledge most CFD efforts have 
been somewhat isolated, so it has been difficult to make an assessment of the state-of-the-art as 
a whole. In order to more broadly assess the current capabilities of different classes of computa- 
tional methodologies, the time appeared to be ripe to hold a workshop for which many participants 
would compute the same flow control test cases. To this end, the CFD Validation of Synthetic 
Jets and Turbulent Separation Control (CFDVAL2004) workshop was held in Williamsburg, Vir- 
ginia in March 2004. 14 Three different test cases, all of which were carried out experimentally at 
NASA Langley Research Center, exercised various aspects related to the flow physics of separation 
control. 

This paper is a summary report of the overall CFD results from the CFDVAL2004 workshop. 
It does not directly address technical issues associated with the experiments themselves, except 
as they relate to the CFD modeling of the flows. Technical papers describing the experiments are 
also available. 15-17 This paper is organized as follows. First, a brief description of the workshop 
is given, including an overall description of the three test cases. Then, a summary of results from 
each of the three test cases is presented. Finally, the paper closes with summary and conclusions. 

II. The CFDVAL2004 Workshop 

The CFDVAL2004 workshop was structured similarly to the series of Refined Turbulence Mod- 
eling workshops sponsored by the European Research Community on Flow, Turbulence, and Com- 
bustion (ERCOFTAC). The 10th such workshop was held in October 2002 in Poitiers, France. 
Although not formally a part of the European series, CFDVAL2004 was held in association with 
the ERCOFTAC organization, and representative members from that organization contributed to 
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the workshop. CFDVAL2004 was also associated with the Air Force Office of Scientific Research 
(AFOSR), the International Association of Hydraulic Engineering and Research (IAHR), QNET- 
CFD, and the National Institute of Aerospace (NLA). 

There were 75 attendees at the workshop. Seven countries were represented, including the 
United States, France, Italy, Germany, Japan, United Kingdom, and Switzerland. Most of the 
workshop participants came from universities, but several companies and public sector research 
laboratories were represented as well. 

The purpose of CFDVAL2004 was to bring together an international group of computational 
fluid dynamics practitioners to assess the current capabilities of different classes of turbulent flow 
solution methodologies to predict flow fields induced by synthetic jets and separation control ge- 
ometries. The experimental data for the workshop were taken specifically with CFD validation in 
mind. Therefore, a great deal of effort was expended to attempt to document most of the relevant 
boundary conditions. The 3-day workshop gave time for each participant to discuss their method 
and results. Each test case was concluded by a summary talk by the organizers, for which all of 
the contributed results were compiled. A panel discussion, whose focus was quality and trust for 
synthetic jet CFD and experiment, was also included in the workshop. 

To encourage broad participation and to determine the general state-of-the-art, the decision was 
made early not to dictate particular boundary conditions, grids, or method of solution. For example, 
although the experimental velocity was given as a function of time near the centers of the jet exits 
for cases 1 and 2, the participants were given the choice as to the precise boundary condition 
used to try to match it. While this strategy to broaden workshop participation was successful - 
solution methods ranged from reduced-order models through Reynolds-averaged Navier-Stokes 
(RANS), to large-eddy simulation (LES), and direct numerical simulation (DNS) - it also had an 
unavoidable downside. Differences in grids and boundary conditions represented an additional 
source of uncertainty when attempting to compare CFD results with each other. 

The three test cases were chosen to represent different aspects of flow control physics. In 
test case 1 (synthetic jet into quiescent air), flow passed in and out of a slot (1.27 mm wide by 
35.56 mm long), which was located on the floor of an enclosed box 0.61 m per side. The jet was 
driven by a side-mounted circular piezo-electric diaphragm inside the cavity chamber beneath the 
floor. The frequency was approximately 445 Hz, and the maximum velocity out of the slot was 
approximately 25-30 m/s. This case was considered to be nominally two-dimensional at the center 
plane of the slot. For case 1, the participants were asked to supply long-time-average jet width and 
velocity profiles at several locations. Also, phase-averaged velocities and turbulence quantities 
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were requested at 8 different phases at several locations. Finally, time-history values of velocity 
at three points in the flow, as well as specific line contour plots of phase-averaged velocities were 
also requested. 

In test case 2 (synthetic jet in a crossflow), flow passed in and out of a circular orifice 6.35 
mm in diameter. The orifice was located on the floor of a wind tunnel splitter plate with a turbu- 
lent boundary layer at M=0.1 and approximate boundary layer thickness of 21 mm. The jet was 
driven electro-mechanically by a bottom-mounted square-shaped rigid piston mounted on an elas- 
tic membrane inside the cavity chamber beneath the splitter plate. The cavity was approximately 
1.7 mm deep and the piston moved approximately ±0.77 mm. The frequency was 150 Hz, and 
the maximum velocity out of the orifice was approximately 1.3 (Too. For case 2, both long-time- 
average jet width and velocities as well as phase-averaged velocities and turbulence quantities (9 
phases) were requested at several locations. Furthermore, participants were asked to supply time- 
history values of velocity at three points in the flow, as well as specific line contour plots of both 
long-time-average and phase-averaged u-velocity. 

In test case 3 (flow over a hump model), turbulent flow at M=0.1 passed over a hump of chord 
420 mm mounted on the floor of a wind tunnel splitter plate. The hump used end-plates at both 
sides (the model had a span of 584.2 mm between the plates), and the flow was nominally two- 
dimensional at its center plane. The hump had a slot near 65% chord, near where separation 
naturally occurred. This case had two mandatory conditions: no-flow-control (no forced flow 
through the slot), and steady suction (m = 0.01518 kg/s). There was also one optional condition 
(the experiment was not completed in time for the workshop) of oscillatory (synthetic) jet control. 
In this case, the jet was driven electro-mechanically by a bottom-mounted rectangular-shaped rigid 
piston mounted on an elastic membrane deep inside the cavity chamber. The frequency was 138.5 
Hz, and the peak velocity out of the slot was approximately 27 m/s. For the case 3 no-flow- 
control and steady-suction conditions, participants were asked to supply Cp and 67/ along the 
hump surface, velocity and turbulence profiles at several locations, and streamline contour plots. 
For the optional synthetic jet condition, only long-time-average Cp along the hump surface was 
requested. 

There was a great deal of CFD data submitted to this workshop. It is clearly not possible 
to show most of the results here. Summary plots as well as several specific detailed plots are 
shown that we believe to be representative of the solutions as a whole, or which serve to illustrate 
specific points or differences. Additional detailed plots of the CFD results are available from 
the CFDVAL2004 website. 14 It should also be noted that workshop participants were allowed to 
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correct, resubmit, or withdraw their submissions during the month following the workshop; both 
this paper summary and the plots on the website contain updated information as of May 2004. 


III. Case 1: Synthetic Jet into Quiescent Air 

A general tabulation of the submissions for Case 1 is given in Table 1 (additional details can 
be found from the CFDVAL2004 website 14 ). There were 8 contributors who ran 25 separate cases. 
There was one LES submission, one reduced-order model submission, one 2-D blended RANS- 
LES submission, and several laminar Navier-Stokes submissions (the 3-D laminar runs perhaps 
can also be characterized as under-resolved DNS); the others were unsteady RANS (URANS). 
Most of the runs were computed in 2-D, but there were also a few 3-D computations. Note that 
none of the 3-D computations modeled the actual shape of the the cavity, including the circular 
diaphragm; instead, they were all computed using periodicity in the direction aligned with the 
slot’s long axis. Six of the contributors modeled the cavity (or some approximation of the cavity) 
and two did not model any cavity. Of those who modeled the cavity, ONERA-flu3m and UKY- 
ghost applied a time-varying velocity on the side of the cavity where the diaphragm was located, 
derived from diaphragm-center displacement data. WASHU-wind and NASA-tlns3d also applied 
a similar boundary condition there, but based it on best matching of the data at the slot exit. GWU- 
vicar3d used a different cavity shape with time-varying velocity boundary condition applied at the 
bottom wall. NCAT-quasld modeled the motion of the actuator with a quasi- 1-D model. 

All of the results were calculated with structured grids. Only two contributors (UKY-ghost 
and NASA-tlns3d) ran with the same grid, “2-D Structured Grid #1”. This grid, along with other 
structured and unstructured grids, were made available on the website several months prior to the 
workshop. Descriptions of these grids can be found at the website. 14 

In this section, we compare results for case 1 with particle image velocimetry (PIV) and hot 
wire data. It should be noted that other experimental data were later taken using a different piezo- 
electric diaphragm (these devices are prone to occasional failure). This later data also include 
laser doppler velocimetry (LDV) measurements in addition to PIV and hot wire, but they were 
also at slightly different conditions (higher maximum outflow velocity) than the data used for the 
workshop. The later data are not shown here, but are reported in Yao et al. 15 A sketch showing the 
locations where most of the results from the computations were supplied is given in Fig. 1, and the 
definition of jet width used for this case is given in Fig. 2. 

Time histories of v- velocities near the center of the slot exit are shown in Fig. 3. There is 
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clearly a significant variation between the two measurement techniques (PIV and hotwire), so it 
is difficult to make a confident assessment of the CFD results. Both experimental measurements 
indicated approximately a 200° phase difference between peak and trough, while many of the CFD 
methods yielded 180°. Some of the differences in particular phase-averaged results between CFD 
and experiment could be a result of not matching the same phase. 

Fig. 4 shows long-time-average v-velocity profiles along the centerline, x=0. Again, a sig- 
nificant difference can be seen between the two measurement techniques, especially very near 
the slot for y less than approximately 2 mm. There were also significant variations between the 
CFD results. Most of them showed similar trends to the experiment, although NCAT-quasld and 
WARWICK-neat indicated a more rapid drop in velocity for y > 10 mm compared to the other 
results. 

Long-time-average results along two y=constant lines (y=0.1 mm and y=4 mm, respectively) 
are shown in Figs. 5 and 6. Just above the slot exit, at y=0.1 mm, all results except POIT-satume 
greatly underpredicted the average jet magnitude. However, beginning near 1-2 mm and above, all 
the CFD results agreed very well with the PIV experiment. An example of this better agreement 
is shown at y=4 mm in Fig. 6. There was relatively little spread among the CFD results at this 
location. Although not shown, the variation between the CFD results increased at higher y-stations. 

The long-time-average jet width is shown in Fig. 7. Most CFD results tended to somewhat 
overpredict the width. The ONERA-flu3m laminar result was similar to the other methods near 
to the wall, but then predicted too much spreading of the jet past y=4 mm. It was more than 
twice as wide as the experimental results near y=8 mm. Although not shown, the CFD results had 
significantly greater variation for y > 8 mm. 

For the sake of brevity, only one representative phase-averaged result is shown here. Fig. 8 
shows v-velocity profiles at y=4 mm, at phase=135°. Compared with the long-time-average results 
shown earlier, this typical phase-averaged result illustrates that there was a great deal more CFD 
variation at this level of detail. Also, the comparison was not as good with experimental data. 

Several participants investigated the effect of time step and grid size on the solution. Although 
not shown, time step had very little effect for any of the solutions. The effect of grid was relatively 
small near the wall, but it could be larger away from the wall. This is illustrated in Fig. 9, which 
shows long-time-average v-velocity along the centerline. UKY-ghost and POIT-saturne showed 
only minor variation, but the 3-D GWU-vicar3d results were affected by grid size and its spanwise 
extent. NASA-tlns3d fine and medium grids gave nearly identical results, but the coarse grid was 
different. 
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The turbulence model’s effect on the CFD results was found to be fairly significant for this case. 
Results are shown in Fig. 10. For POIT-saturne, the k-e model compared better with experimental 
data than the Reynolds stress model (RSM) for this quantity. For WARWICK-neat, three different 
turbulence models were essentially the same for y < 4 mm, but then they behaved very differently 
beyond that. None of these agreed particularly well with experimental data. For WASHU-wind, 
Spalart-Allmaras (SA) and Menter’s shear stress transport (SST) were fairly similar, especially for 
y < 7 mm, but the 2-D blended SST-LES method generally overpredicted the average velocity 
everywhere. For NASA-tlns3d, SA and SST showed similar variation to that seen by WASHU- 
wind: the two models were close near the wall, then SST predicted somewhat larger velocities 
than SA further from the wall. 

Although not shown, the computed turbulence quantities (e.g., u'v') did not compare well at 
all with experimental data, and there was very large variation exhibited among CFD results. Also, 
there is some question whether the measured data is entirely a result of turbulence, or a result 
of unsteady “flapping” of the time-dependent jet. It was generally felt by the participants of the 
workshop that this case 1 flow is probably laminar in the beginning, then transitional further into the 
flow field. Therefore the use of a standard RANS turbulence model may not be entirely warranted. 

As mentioned earlier, in spite of attempting to match the v-velocity phase at the slot outflow 
(see Fig. 3), many of the participants had difficulty because the experimental variation was not 
sinusoidal (180° peak-to-peak). As a result, at certain phases they found that the computed jet 
could be at a different position than in the experiment. This is illustrated by example flowfield 
contours at phase=135°, shown in Fig. 11. Similar to most other results (not shown), the NASA- 
tlns3d result at this phase had the extent of its peak jet velocity lower (near 4 mm) than seen in 
the experimental data (near 7 mm). It is also possible that three-dimensional effects due to ring 
vortices formed from the slot ends had an influence on the flow field in the experiment, especially 
at distances far from the wall. 

For this case, the 2-D laminar ONERA-flu3m result predicted too much jet spreading. How- 
ever, the 2-D laminar NCAT-quasld, 3-D laminar, and 3-D LES results were consistent with other 
CFD methods in this regard. Different turbulence models could have significant effects, especially 
further into the flow field: the SA, SST, and POIT-saturne -ke models seemed to be of similar qual- 
ity, and generally agreed best with the data among the URANS models. Surprisingly, there did 
not appear to be a distinct advantage to modeling the cavity for this case, as opposed to specifying 
a jet boundary condition at the surface. Regarding grid effect, 2-D grid levels of about 60,000 
points (including cavity) were found to be sufficient to capture the flow physics, whereas 16,000 
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points yielded noticeable differences in the solution. Among 3-D results, changing from 460,000 
to 700,000 points (or increasing the spanwise periodic extent) could significantly alter the solution 
far into the flow field. There was very little effect of changing time step, even with as low as 72 
steps per period. 


IV. Case 2: Synthetic Jet in a Crossflow 

A general tabulation of the submissions for Case 2 is given in Table 2 (additional details can 
be found from the CFDVAL2004 website 14 ). There were 5 contributors who ran 10 separate cases. 
There was one LES submission; the others used URANS. All runs were (necessarily) 3-D. All 
methods were second order in space and time. All submissions modeled the cavity, except for 
CIRA-zen which specified a time-varying profile at the orifice exit. However, both USTO-rans 
and NASA-fun3d altered the cavity shape. In the latter case, the bottom wall was made flat while 
keeping the volume the same as the original (stepped) geometry. Two of the contributors modeled 
the full plane, while the others modeled a half plane with symmetry imposed at the center plane. 
Each contributor used a different grid: four used structured grids and one used unstructured. The 
“Structured Grid #1” and “Unstructured Grid #1” were made available on the website several 
months prior to the workshop. Descriptions of these grids can be found at the website. 14 The 
four contributors who modeled the cavity applied a time-dependent velocity specification at the 
cavity bottom wall: V = Vcos(2n ft). Each adjusted V to achieve what was considered to be a 
reasonable match of the experimental w-velocity variation at the orifice exit. 

In this paper, case 2 results are quantitatively compared only to LDV data. A sketch showing 
some of the LDV data measurement locations is given in Fig. 12. PIV data was also acquired at 
several planes. The reader interested in quantitatively comparing experimental results using the 
two different techniques is referred to Schaeffler and Jenkins. 16 

Time histories of velocities near the center of the orifice exit are shown in Fig. 13. The u- 
velocity component was captured well by three of the codes. CIRA-zen showed a higher value 
and NASA-fun3d showed a lower value than experimental results during the expulsion part of the 
cycle. All codes predicted a zero or near- zero cross-stream v-component of velocity, whereas the 
experimental data indicated a very large v-component during expulsion. The cause of this high 
experimental value is not known. Finally, all codes captured the w-velocity reasonably well (with 
some variation), although no one replicated the “dip” near phase=160°. The cause for this dip 
in the experimental result is also not known. The cavity under the orifice was extremely shallow 
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(nominally 1.7 mm deep neutral position with tunnel on), and the piston moved ±0.77 mm up and 
down from this position. The volume in the cavity more than doubled from roughly 10,000 to 
26,000 cubic mm during the cycle. It is questionable whether a simple non-moving-wall transpi- 
ration boundary condition, used by all of the CFD methods, is a very accurate model for this large 
a range of piston motion. 

Fig. 14 shows long-time-average u-velocity profiles on the center plane at three downstream 
stations. Overall, there was fair agreement between CFD and experiment, but there was also a 
significant spread among the CFD results. Fig. 15 shows w-velocity profiles along one spanwise 
line and two streamwise lines. Along the spanwise line, there were no experimental LDV data, and 
there was a significant spread in the CFD results. Along the streamwise lines, the 5 CFD results 
again exhibited significant spread; two of the CFD methods exhibited reasonably good agreement 
with experimental results along the z=10mm, y=0 line (Fig. 15c). 

Two of the participants examined the effect of grid size on the solution. The results can be 
briefly summarized in a representative figure, Fig. 16. In general, the fine (3.9 million cells) and 
coarser (0.49 million cells) grids used in NASA-cfl3d did not exhibit a dramatic influence in the 
solution. However, for NASA-fun3d the fine (0.26 million nodes) and coarser (46,000 nodes) grids 
did have a significant difference. For NASA-fun3d, the finer grid produced a negative streamwise 
velocity during expulsion, but the coarser grid did not. 

Two of the participants investigated the effects of different turbulence models on the solution. 
Again, only one representative figure is shown here for the sake of brevity, Fig. 17. In NASA- 
cfl3d, three different models resulted in minor differences from each other, and none was clearly 
better in comparison with the experimental data. The two models used by USTO-rans also showed 
only minor differences from each other, and they were clearly different in character from both the 
NASA-cfl3d results and the experimental data. In other words, the turbulence models themselves 
did not have as much of an impact as the combined effect of different codes, grids (in number of 
points and full plane vs. half plane), and other solution variations (such as alteration of the cavity 
shape by USTO-rans). 

Only a few phase- averaged results are shown here. Fig. 18 shows u-velocity along a line ID 
downstream at three different phases. All CFD results were in good agreement with each other 
and with experimental results at phase=0° (at this time the orifice is still in its suction phase, so 
there is little to no dynamic flow occurring ID downstream). However, at the two later phases the 
influence of the expulsion part of the cycle is felt at this location. The CFD results captured the 
influence, generally following the trends seen in the experiment, but there was a significant spread 
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in the solutions. 

Phase-averaged turbulence quantities from the CFD solutions generally exhibited similar trends 
to each other, but again with a significant spread. There were usually larger discrepancies between 
CFD and experiment than the discrepancies seen between turbulence profiles. An example is 
shown in Fig. 19. 

Finally, a sample comparison between experiment and CFD over the entire plane ID down- 
stream is shown in Fig. 20. This figure shows u-velocity contours at phase=120°. It is intended as 
a qualitative comparison only. Also, only one CFD result is shown for brevity; other CFD results 
were qualitatively similar. Although details were different, this figure indicates that CFD predicted 
the overall structure of the dynamic flow resulting from the expulsion part of the cycle as it passes 
downstream. 

For this case study, LES vs. URANS did not appear to be a significant factor. In fact, overall 
(in a subjective sense) the two methods that tended to yield the most similar long-time-average and 
phase-averaged mean-flow solutions were ONERA-flu3m (LES) and NASA-cfl3d (URANS). In 
general, no one method, algorithm, or turbulence model stood out as being the best methodology 
for capturing the physics of this flow. Grid levels of about a half a million points or more appeared 
to be sufficient to capture the general character of the unsteady flow sufficiently well, whereas 
grids coarser than this could change the solution significantly. None of the participants submitted 
results examining the effect of time step; however, NASA-cfl3d performed computations that were 
not submitted, varying between 720 and 1440 steps per cycle, and found very little effect. 

V. Case 3: Flow over a Hump Model (Separation Control) 

A general tabulation of the submissions for Case 3 is given in Table 3 (additional details can be 
found from the CFDVAL2004 website 14 ). There were 13 contributors who ran 56 separate cases. 
Methods were mostly RANS/URANS, but there was also one DNS result and several blended 
RANS-LES results. Most of the runs were computed in 2-D, but there were also several 3-D com- 
putations. Some of the 3-D computations were periodic in the spanwise direction, and some were 
half-plane (either flat inviscid side wall or actual side plate geometry modeled viscously). Most of 
the contributors modeled the cavity, but several applied boundary conditions directly on the hump 
surface. Results were computed on both structured and unstructured grids. Many contributors used 
the 2-D grids supplied on the website (for example, “2-D Structured Grid #1”). Descriptions of 
these grids can be found there. 14 
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Case 3 is different from cases 1 and 2 in that most of the results were not time-dependent. 
Only the optional oscillatory-control case was unsteady, and for that case only long-time-average 
Cp’s were requested from the participants (8 of the contributors computed this, running 12 sepa- 
rate cases). However, when there is more extensive time -dependent data available in the future, the 
oscillatory-control experiment may serve as a good test case for time-accurate validation exercises. 
A sketch showing the locations where most of the results from the computations were supplied is 
given in Fig. 21. Computed pressures are compared with experimental data from surface-mounted 
pressure taps, and velocity and turbulence profiles are compared with PIV experimental data. Ex- 
perimental details are reported in Greenblatt et al. 17 

A summary of all of the CFD results for the no-flow-control condition is shown in Fig. 22. 
(Note that separation and reattachment locations are not displayed in the figures for those methods 
that did not report skin friction.) As a whole, most CFD results missed the pressure levels over 
the hump between 0.2 < x/c < 0.6, and also predicted higher pressures than experimental results 
in the separated region upstream of x/c=l. As will be shown below, the missed pressures were 
likely due to blockage effects caused by the side plates in the experiment. The separation location 
(not known precisely for the experiment) was predicted reasonably well by most CFD methods, 
and the reattachment location was predicted (for the most part) significantly downstream of the 
experimental location of x/c = 1.11. Most models of similar type generally behaved similarly. 
Two exceptions were CIRA-zen-sst-kprof-3, which predicted separation later and reattachment 
earlier than other results with the SST model, and CIRA-zen-ke, which predicted reattachment 
further downstream than other results using a k-e model. CIRA-zen also reported unusually low 
skin friction levels for all of its results. 

A summary of all of the CFD results for the steady suction condition is shown in Fig. 23. Again, 
most CFD results overpredicted the pressure levels over the hump between 0.2 < x/c < 0.6. Also, 
the pressures in the separated region did not exhibit the same levels or shape as experimental 
results. The separation location (again not known precisely for the experiment) was predicted to 
be slightly upstream of the experimental location, and the reattachment location was predicted 
to be significantly downstream of the experimental location of x/c = 0.94 by all CFD methods. 
For the suction condition, there was generally a broader range exhibited in the results for a given 
turbulence model than for the no-flow-control condition. Notable results that showed the largest 
differences from other turbulence models of the same type were: CTR-fluent-sa-2, CTR-fluent-sa- 
1, NASA-rans-sst-weno5, and US-fluent-ke. 

The two figures, Figs. 22 and 23, are also instructive from the point of view of CFD uncertainty. 
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The surface pressure coefficient plots illustrate the range of variation of CFD for these cases. This 
variation is caused by use of different grids, codes, turbulence models, and boundary conditions 
(including 2-D vs. 3-D modeling). In particular, it should be noted that the downstream boundary 
condition in internal flow cases such as this can affect the flow to some degree: by setting different 
back pressures, one can achieve different results. Individual differences in this boundary condition 
likely account for some of the CFD variation. 

The tendency for most of the CFD results to miss the pressure levels over the front half of 
the hump was believed to be due to blockage effects in the experiment. Two pieces of evidence 
support this conjecture. First, the 3-D AZ-cobalt runs that accounted for the side plate physical 
shape (modeled as a viscous wall) resulted in improved levels there. Second, the 2-D NASA- 
cfl3d run that used an altered top wall shape to approximately account for the side-plate blockage 
(using area ratio of plate cross-section to tunnel cross section) resulted in excellent agreement of 
pressures over the front part of the hump compared with experimental data. These results are 
shown in Fig. 24. The grid used in the NASA-cfl3d case (2-D Structured #5) is described on the 
CFDVAL2004 website; this grid was created after the workshop. It should be noted from this 
figure that the 3-D AZ-cobalt runs also yielded good results in the separated region, whereas the 
2-D NASA-cfl3d runs did not. This suggests that there may be 3-D separated structures along the 
back end of the side plates that further constrict the flow (additional blockage) in the separated 
region. 

Although not shown, most of the contributed results were in fairly good agreement with ex- 
perimental data (for both velocity profiles and turbulent shear stress profiles) just upstream of the 
slot (especially for the no-flow-control condition). And the velocity profiles inside the separated 
region were also predicted reasonably well, on the whole, for both no-flow-control and suction 
conditions. However, because the CFD results reattached too late, most velocity profiles compared 
poorly with experimental data at stations downstream of reattachment. The two exceptions to this 
were AZ-cobalt-des-l-3d and META-cfd++lns-3d for the no-flow-control condition. Profiles us- 
ing these methods (along with five other typical representative solutions) at x/c = 1.2 are shown 
in Fig. 25. However, all results were poor downstream of reattachment for the suction condition. 
Fourteen representative results are shown at x/c = 1.0 in Fig. 26. 

A possible reason for reattachment being predicted too late is that most of the current models 
and methods predicted turbulent shear stress to be too small in magnitude in the separated region. 
As an illustration of this tendency, Fig. 27 shows both u-velocity and turbulent shear stress for the 
suction condition at the location x/c = 0.8, using several representative solutions. The velocity 
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profiles exhibited a variation in shape, especially near the center of the profile, but the overall agree- 
ment with experimental results was fairly good for the RANS methods. However, the magnitude of 
the turbulent shear stress was seriously underpredicted. (The turbulent shear stress data had a max- 
imum error of as much as 20%, based on a pseudo-empirical uncertainty analysis. Nonetheless, 
the results still indicated a large difference between CFD and experiment.) One of the anomalies 
evident in the submissions was that the US-fluent results consistently yielded significantly larger- 
in-magnitude turbulent shear stress levels, in better agreement with experimental data, even though 
they used many of the same turbulence models as other codes. In any case, the US-fluent solutions 
still predicted reattachment too far downstream, similar to the other submissions. 

The DNS and blended RANS-LES solutions (DES, LNS, FSM) were generally on a par with 
the RANS methods; they, too, predicted reattachment too far downstream. However, DNS, LES, 
and 3-D blended RANS-LES models involve typically finer grids than RANS methods, and require 
very lengthy run times. It is possible that future improvements may result from additional efforts in 
these areas. As shown earlier in Fig. 25, a few of the blended RANS-LES models appeared to do a 
better job predicting velocity profiles near and downstream of reattachment than the RANS meth- 
ods for the no-flow-control condition, but they did not improve results for the suction condition. 
Also, it is worth noting that AZ-cobalt-des predicted resolved turbulent shear stress levels down- 
stream of reattachment (for both no-flow-control and suction conditions) that were much higher in 
magnitude than most RANS results. 

Many of the contributors investigated the effect of grid on the solution, for the no-flow-control 
and suction conditions. For 2-D, little influence was noted (all used a minimum of about 50,000 
grid points). The influence of turbulence model has been shown earlier in Figs. 22 and 23. Gener- 
ally, there were some differences between the models as expected, but all performed poorly when 
it came to predicting reattachment location. No one model stood out as performing consistently 
better for this case. A few contributors ran without modeling the cavity. In general this did not 
appear to either help or hurt, even for the suction case. CIRA-zen found that applying oblique 
suction improved results over applying normal suction. 

A summary of all of the CFD results for the (optional) oscillatory control case is shown in 
Fig. 28. There was variation among the CFD results, but in general they all exhibited similar 
trends. The results whose separated region pressure recovery occurred furthest upstream (in better 
agreement with experimental data) were UAZ-cfl3d, BOEING- overflow, and UMD-rans. NASA- 
fun2d investigated the effect of time step on this solution, and found very little effect. 
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VI. Summary and Conclusions 


Time-dependent flows involving unsteady flow control are difficult not only to compute with 
CFD, but also to measure experimentally. The CFDVAL2004 workshop in many ways ended up 
being as much of a workshop on experimental issues as on CFD issues. From the two cases involv- 
ing synthetic jets (case 1 and case 2), we learned that there are still some inconsistencies between 
different measurement techniques for certain aspects of these types of time-dependent flows. These 
differences were openly reported at the workshop, in order to foster honest and open dialog and 
to help encourage collaborative efforts to work toward improving future results. Most importantly, 
the workshop proved to be an ideal setting for providing a quality technical interchange between 
CFD and experimental scientists, where each could learn about the capabilities and limitations 
inherent in the processes and tools of the other. 

This paper focused on providing a summary of the CFD results from the workshop. The work- 
shop itself was designed to encourage wide participation and to establish a baseline for the state- 
of-the-art in the area of synthetic jet validation and turbulent separation control. Specific grids, 
boundary conditions, and methods of solution were not dictated to the participants. As a result, 
CFD variation was fairly large. However, trends could still be noted, including overall capabilities 
and shortcomings. Many of the workshop participants suggested that, now that we have estab- 
lished a baseline for these cases, future CFD efforts should attempt to minimize known sources of 
CFD variation, especially regarding boundary conditions. It was also suggested that future exper- 
iments focus additional effort documenting the time-dependent boundary conditions, especially at 
and near the exit plane of the jet/suction slot or orifice. 

The bottom line from the CFD results can be summed up as follows: no one CFD technique ex- 
celled above others, and there was wide variation (especially for time-dependent results), and only 
qualitative agreement with experimental data. In other words, the “state-of-the-art” CFD methods 
of today are not fully adequate to consistently and accurately predict these types of flows. What 
is today’s “state-of-the-art”? Most of the contributors used URANS for the time-dependent cases 
and RANS to compute the steady cases. A wide variety of turbulence models were employed. 
Most CFD methods were second order in space and time. Some contributors used higher order, 
but there did not appear to be any obvious benefits from doing so for these cases. Most contrib- 
utors utilized transpiration boundary conditions on a fixed grid to simulate moving walls. The 
few blended RANS-LES, LES, and DNS solutions showed merit, and were of similar quality to 
the RANS and URANS solutions. We believe that with today’s increased computer power and 
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improved algorithms, these advanced methods are currently coming of age and represent the way 
much of CFD will be performed in the future. However, as yet they showed no clear benefits over 
RANS/URANS in the sense of providing consistently better results for the workshop cases. 

Case 1 (synthetic jet into quiescent air) was a difficult experiment to simulate. The flow field 
was probably mostly laminar or transitional, so it was unclear how best to simulate it. The piezo- 
electric driver used in the experiment was difficult to model using CFD because it moved like 
a drum (and possibly with multiple shape modes), and not like a rigid piston. The experiment 
indicated a deviation from periodicity in the velocity near the jet exit that was not completely 
simulated or captured by the CFD results, and caused specific phase results to be mis-aligned. For 
example, at phase=135° in the experiment, the jet extended higher into the field than at the same 
phase as defined by most of the CFD runs. Long-time-average results showed deviations from 
experimental data very near to the slot as well as far from the slot. In the intermediate range (4-8 
mm from the slot), the results were generally better. CFD yielded significantly larger variations 
in the phase-averaged quantities. Except for a few cases, average jet width was usually somewhat 
overpredicted by CFD in the near field. In the URANS simulations, turbulence models were found 
to have a significant effect for this case. In general, the k-e, SA, and SST models performed better 
compared to experimental data than other models tested. 

Case 2 (synthetic jet in a crossflow) was the least computed of the three cases, probably owing 
to its being both time-dependent and (necessarily) 3-D. Similar to case 1, the time-dependent 
experimental velocities measured near the orifice exit exhibited anomalies not captured or modeled 
by CFD. In particular, the experiment exhibited a dip in w-velocity on the downstroke and a very 
large spanwise velocity component during the expulsion part of the cycle that none of the CFD 
methods accounted for. In spite of this, reasonably good qualitative results were obtained compared 
to experimental results, but there were significant CFD variations in both the long-time-average 
and the phase-averaged results. No one method, algorithm, or turbulence model stood out as being 
the best methodology. Even LES vs. URANS did not appear to be a significant factor. Generally, 
different turbulence models for URANS did not have as much of an impact on this case as different 
grids, codes, and other solution variations. 

Case 3 (flow over a hump model) had the greatest number of workshop participants, since the 
required cases were nominally steady and 2-D. However, it was discovered after the workshop that 
the side plates used in the tunnel caused blockage that, if not modeled, resulted in relatively mi- 
nor (but noticeable) overprediction of the pressures over most of the hump. Also, based on some 
3-D results that modeled the side plates viscously and obtained good pressure comparisons in the 
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separated region, it is likely that flow structures near the back of the plates constrict the flow even 
further in that region. In other words, in spite of this flow appearing to be relatively simple, com- 
puting the wall pressures accurately can require full 3-D modeling or else some sort of blockage 
corrections. Overall, case 3 was not a severe test from the point of view of predicting separa- 
tion. But CFD results were deficient in another important regard: they consistently predicted the 
reattachment location to be significantly further downstream than the location documented in the 
experiment. This same behavior occurred regardless of turbulence model or method; even a DNS 
computation predicted too long of a separation bubble. Inside the bubble itself, most computations 
predicted velocity profiles in reasonably good agreement with experimental data, but turbulent 
shear stresses were underpredicted in magnitude. This underprediction is consistent with delayed 
reattachment, because it indicates reduced turbulent mixing inside the separated region compared 
with experimental data. 

In conclusion, the CFDVAL2004 workshop has established a benchmark for three different 
flows involving synthetic jets and turbulent separation control. Although the current state-of-the- 
art CFD methods are deficient in being able to consistently and accurately predict these flows, two 
areas have been identified as key to improvement. 

• First, for synthetic jets, it is important to employ consistent boundary conditions when mul- 
tiple CFD methods are compared. Only then will we be able to isolate and work to correct 
deficiencies in the CFD models and methods. Implicit in this need is the requirement that 
experiments be used to document extremely detailed and accurate time-dependent flow field 
variables at and near the slot/orifice exits. 

• Second, for the hump model case, turbulence models (for RANS) or other methods such 
as LES or blended RANS-LES need to be improved/developed/calibrated to increase the 
mixing in the separated region and bring about earlier reattachment and flow recovery. 
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Table 1. Case 1 Summary of Submissions 


Label 

Organization 

Authors 

Method 

ONERA-flu3m 

ONERA 

Mary 

LES, Lam. N-S, URANS (SA) 

UKY-ghost 

U. Kentucky 

Huang, Chen, Katam, 
Parimi, LeBeau, Huang 

URANS (SST) 

GWU-vicar3d 

GWU 

Rupesh, Ravi, Mittal, 
Raju, Gallas, Cattafesta 

Lam. N-S 

NCAT-quasld 

NC A&T State U. & 
NASA LaRC 

Yamaleev, Carpenter, 
Vatsa 

Reduced-order mod. 
in slot+4th ord. Lam. N-S 

POIT-saturne 

U. Poitiers 

Carpy, Manceau 

URANS (k-e, DRSM) 

WARWICK-neat 

U. Warwick & U. Wales 

Preece, Tucker, Liu 

URANS (k-e, nonlin. k-e, EASM) 

WASHU-wind 

Washington U. 

Cui, Agarwal 

URANS (SA, SST, SST-LES) 

NASA-tlns3d 

NASA LaRC 

Vatsa, Turkel 

URANS (SA, SST) 


Table 2. Case 2 Summary of Submissions 


Label 

Organization 

Authors 

Method 

NASA-cfl3d 

NASA LaRC 

Rumsey 

URANS (SA, SST, EASM) 

USTO-rans 

USTO & ETH 

Azzi, Lakehal 

URANS (TLV, EASM) 

ONERA-flu3m 

ONERA 

Dandois, Gamier, Sagaut 

LES 

CIRA-zen 

CIRA 

Marongiu 

URANS (k-e) 

NASA-fun3d 

NASA LaRC 

Atkins 

URANS (SA) 
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Table 3. Case 3 Summary of Submissions 


Label 

Organization 

Authors 

Method 

AFRL-fdl3di 

OAI & AFRL 

Morgan, Rizzetta, Visbal 

RANS fk-e) 

AZ-cobalt 

Arizona State & 
Cobalt Solutions 

Krishnan, Squires, 
Forsythe 

DES, RANS (SA, SST) 

BOEING-overflow 

Boeing 

Shmilovich, Yadlin, Clark 

RANS (SA, SST) 

META-cfd++ 

Metacomp Tech 

Shariff, Batten 

RANS (cubic k-e), LNS 

CTR-fluent 

CIRA-zen 

CIRA & CTR 

Marongiu, Iaccarino, 
Catalano, Amato 

RANS (SA) 
RANS (SST, k-e) 

NASA-cfl3d 

NASA LaRC 

Rumsey 

RANS (SA, SST, EASM) 

NASA-rans 

NASA LaRC 

Balakumar 

RANS (SST) 

NASA-fun2d 

NASA LaRC 

Viken 

RANS (SA) 

UAZ-cfl3d 

U. Arizona 

Israel, Fasel 

FSM w/lin. EASM 

UAZ-DNS 

U. Arizona 

Postl, Wernz, Fasel 

DNS 

UK-ghost 

U. Kentucky 

Katam, Chen, Huang, 
Parimi, LeBeau, Huang 

RANS (SST) 

UMD-rans 

U. Maryland 

Duraisamy, Baeder 

RANS (SA, SST) 

US-fluent 

Utah State 

Spall, Phillips, Alley 

RANS (k-o, k-e, SA 
SST, v2f) 
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Figure 1. Sketch of the case 1 comparison locations. 
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Figure 2. Sketch showing definition of jet width for case 1. 
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Figure 3. Case 1 time histories of v-velocity at x=0, y=0.1mm. 
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Figure 4. Case 1 average v-velocity profiles at x=0. 


23 of 45 



10 


ONERA-flu3m-les-3d 

ONERA-flu3m-lam 

ONERA-flu3m-sa 

GWU-vicar3d-3d(fine) 

NCAT-quasId+rans 

□ exp, PIV 



J I I I I L 

0 2 

x, mm 


-L 

4 


J 

6 



Figure 5. Case 1 average v-velocity profiles at y=0.1mm. 
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Figure 6. Case 1 average v- velocity profiles at y=4mm. 
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avg jet width, mm 


Figure 7. Case 1 jet width based on average v-velocity. 
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Figure 8. Case 1 phase-averaged v-velocity profiles at y=4mm, phase=135°. 
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Figure 9. Grid effect on case 1 average v-velocity profiles at x=0. 
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Figure 10. Turbulence model effect on case 1 average v-velocity profiles at x=0. 
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Figure 11. Case 1 contours of v-velocity, phase=135°; (a) PIV experiment, (b), NASA-tlns3d-sa(fine). 
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Figure 13. Case 2 time histories at x=50.63mm, y=0, z=0.4mm; (a) u-velocity, (b) v- 
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ity, (c) w-velocity. 
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Figure 14. Case 2 average u-velocity at y=0; (a) ID downstream, (b) 2D downstream, (c) 8D downstream. 
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Figure 15. Case 2 average w-velocity; (a) z=10mm ID downstream, (b) z=0.4mm and y=0, (c) z=10mm and y=0. 
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Figure 16. Grid effect on case 2 u-velocity time histories at x=50.63mm, y=0, z=0.4mm. 



Figure 17. Turbulence model effect on case 2 u-velocity time histories at x=63.5mm, y=0, z=10mm. 
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Figure 18. Case 2 phase-averaged u-velocity at y=0 and ID downstream; (a) phase=0°, (b) phase=120°, (c) 
phase=240°. 
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Figure 19. Turbulence model effect on case 2 phase-averaged u'w' turbulent stress at x=44.45mm, y=0, 
phase=0°. 
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Figure 22. Summary of case 3 results compared to experimental data for no-flow-control condition; (a) Cp, (b) 
separation location, (c) reattachment location. 
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Figure 23. Summary of case 3 results compared to experimental data for suction condition; (a) Cp, (b) separa- 
tion location, (c) reattachment location. 
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Figure 24. Surface pressure coefficients for case 3 results that accounted for tunnel blockage. 


No-flow-control, x/c=1 .2 



Figure 25. Sample of case 3 no-flow-control condition u-velocity profiles at x/c=1.2 (downstream of experimental 
reattachment). 
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Figure 26. Sample of case 3 suction condition u-velocity profiles at x/c=1.0 (downstream of experimental reat- 
tachment). 


43 of 45 


Suction, xl c=0.8 



Suction, x/c=0.8 



Figure 27. Sample of case 3 suction condition results at x/c=0.8 (inside separation bubble); (a) u-velocity, (b) 
turbulent shear stress. 
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Figure 28. Summary of case 3 long-time-average Cp results for oscillatory control condition. 
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